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Abstract 

A first measurement is reported of the longitudinal proton structure function Fl{x, Q"^) at 
the ep collider HERA. It is based on inclusive deep inelastic e^p scattering cross section 
measurements with a positron beam energy of 27.5 GeV and proton beam energies of 920, 
575 and 460 GeV. Employing the energy dependence of the cross section, Fi is measured 
in a range of squared four-momentum transfers 12 < < 90 GeV^ and low Bjorken x 
0.00024 < X < 0.0036. The Fl values agree with higher order QCD calculations based on 
parton densities obtained using cross section data previously measured at HERA. 
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1 Introduction 



This letter presents the first measurement of the longitudinal structure function Fl{x, Q"^) of the 
proton at low Bjorken x. The inclusive deep inelastic ep scattering (DIS) cross section at low 
(5^, written in reduced form as 

is determined by two structure functions, F2 and F^. Here, = —q^ is the negative four- 
momentum squared transferred between the electron^ and the proton, and x = Q'^/2qP denotes 
the Bjorken variable, where P is the four-momentum of the proton. The two variables are 
related through the inelasticity of the scattering process, y = Q'^/sx, where s = AEcEp is the 
centre-of-mass energy squared determined from the electron and proton beam energies, and 
Ep. In equation 1, a denotes the fine structure constant and — 1 -\- [1 — yY. 

The two proton structure functions F^ and F2 are of complementary nature. They are related 
to the 7*p interaction cross sections of longitudinally and transversely polarised virtual photons, 
ul and gt, according to oc ctl and F2 oc ((7^-1- ut). Therefore the relation Q < Fl < F2 
holds. In the Quark Parton Model (QPM), F2 is the sum of the quark and anti-quark x distribu- 
tions, weighted by the square of the electric quark charges, whereas the value of Fl is zero [1]. 
In Quantum Chromodynamics (QCD), the longitudinal structure function differs from zero, re- 
ceiving contributions from quarks and from gluons [2]. At low x and in the region of deep 
inelastic scattering the gluon contribution greatly exceeds the quark contribution. Therefore 
Fl is a direct measure of the gluon distribution to a very good approximation. The gluon dis- 
tribution is also constrained by the scaling violations of F2(,x, Q^) as described by the DGLAP 
QCD evolution equations [3]. An independent measurement of Fl at HERA, and its compar- 
ison with predictions derived from the gluon distribution extracted from the Q"^ evolution of 
F2{x, Q"^), thus represents a crucial test on the validity of perturbative QCD at low Bjorken x. 

The longitudinal structure function, or the equivalent cross section ratio R — aL/<JT — 
FlI{F2 — Fl), was measured previously in fixed target experiments [4] and found to be small 
at large x > 0.2, confirming the QPM prediction in the region of DIS. 

From experimental determinations by HI [5-7], which used assumptions on the behaviour 

of F2 in extracting Fl , and from theoretical analyses of the inclusive DIS cross section data 
[8, 9], the longitudinal structure function at low x is expected to be significantly larger than 
zero. This prediction relies on perturbative QCD calculations of Fl to next-to-leading order 
(NLO) [10] andNNLO [11]. 

The measurement of Fl requires several sets of DIS cross sections at fixed x and but at 
different y. This was achieved at HERA by variations of the proton beam energy whilst keeping 
the lepton beam energy fixed. The sensitivity to Fl is largest at high y as its contribution to is 
proportional to y^. At low Q"^, high y values correspond to low values of the scattered electron 
energy. Small energy depositions can also be caused by hadronic final state particles leading to 

'The term electron is used here to denote both electrons and positrons unless the charge state is specified 
explicitely. The data analysed are from positron-proton scattering, except for some measurements of background 
properties which additionally include electron-proton scattering data. 
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fake electron signals. These are dominantly due to photoproduction processes at Q"^ ~ 0. The 
large size of this background makes the measurement of Fi,{x^ Q"^) particularly challenging. 

The present measurement of Fl{x, Q^) is based on data collected with the HI detector in 
e+p collisions from January to June 2007 with a positron beam energy of 27.5 GeV. Three pro- 
ton beam energies were used, the largest, nominal energy of 920 GeV, the smallest energy of 
460 GeV and an intermediate energy of 575 GeV, chosen for an approximately equal span be- 
tween the three resulting cross section measurements in /Y^ (see equation 1). The integrated 
luminosities collected with HI are 21.6 pb~^, 12.4 pb~^ and 6.2 pb~^, respectively. This letter 
presents first results on Fl in an intermediate range of Q^, between 12 and 90 GeV^. 



2 Data Analysis 
2.1 HI Detector 

The HI detector [12] was built and upgraded for the accurate measurement of ep interactions 
at HERA. The detector components most relevant to this measurement are the central jet drift 
chamber (CJC), the central inner proportional chamber (CIP), the backward lead-scintillator 
calorimeter (SpaCal) and the liquid argon calorimeter (LAr). The CJC measures transverse mo- 
menta of tracks with an accuracy of Spt/pj — 0.005/GeV Complementary tracking information 
is obtained from the backward silicon tracker (BST), which is positioned around the beam pipe, 
and from the z drift chamber COZ, which is located in between the two cylinders of the CJC. 
The CIP provides trigger information on central tracks [13]. The SpaCal [14] has an energy 
resolution of 5E/E ~ 0.07/ \/E/GcV for electromagnetic energy depositions and is comple- 
mented by a hadronic section. It also provides a trigger down to 2 GeV energy. The LAr allows 
the hadronic final state to be reconstructed with an energy resolution of about 0.50/-^£^/GeV. 

Photoproduction events can be tagged with an electron calorimeter placed at 2; = —6m 
downstream in the electron beam direction, which defines the negative z axis and thus the 
backward direction. The luminosity is determined from the Bethe-Heitler scattering process, 
which is measured using a photon calorimeter at 2; = —103 m. 



2.2 Kinematic Reconstruction and Event Selection 

The DIS kinematics at large y are most accurately reconstructed using the polar angle, 9^, and 
the energy, E'^, of the scattered electron according to 

y = 1 - f sin^(^./2) , = , (2) 

where x = Q'^/sy. The event signature of this analysis comprises an electron scattered back- 
wards and a well reconstructed event vertex. The scattered electron energy is measured in the 
backward calorimeter SpaCal. The polar angle is determined by the positions of the interaction 
vertex and the electron cluster in the SpaCal. 
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Energy E'^ of scattered electron candidate 


> 


3.4 GeV 


Transverse size Riog of candidate cluster 


< 


5 cm 


Hadronic energy fraction behind the cluster 


< 


15%0f Ee' 


Transverse distance between cluster and linked track 


< 


6 cm 


E-p, 


> 


35 GeV 


z position of interaction vertex 




J 1 < 35 cm 



Table 1: Criteria applied to select DIS events at high inelasticity y. 



In order to trigger on low energy depositions with a threshold of 2 GeV, a dedicated trigger 
was developed based on the SpaCal cell energy depositions. At small energies the SpaCal 
trigger is complemented by the CIP track trigger which reduces the trigger rate to an acceptable 
level. The efficiency of this high y trigger is constant at around 98 % above 3 GeV, as monitored 
with independent triggers. At energies larger than 7 GeV no track condition is used in the trigger 
and the efficiency, up to highest energies, exceeds 99 %. 

The event selection is based on the identification of the scattered electron as a localised 
energy deposition (cluster) of more than 3.4 GeV in the SpaCal. Hadrons, dominantly from 
photoproduction but also from DIS, may also lead to such energy depositions. This fake electron 
background is reduced by the requirement of a small transverse size of the cluster, Riog, which 
is estimated using a logarithmic energy weighted cluster radius. The background is further 
reduced by the requirement that the energy behind the cluster, measured in the hadronic part 
of the SpaCal, may not exceed a certain fraction of E'^. For lower energies the selected cluster 
must be linked to a track. If the highest energy cluster fails to fulfill the selection criteria, the 
next to highest energy cluster passing the selection criteria is considered. Alternatively ordering 
the SpaCal clusters according to the scattering angle or transverse momentum gives consistent 
cross section results. 

An additional suppression of photoproduction background is achieved by requiring longitu- 
dinal energy-momentum conservation using the variable 

E-p, = Ei{Ei-p,,i) + E'^{1 - cos9e), (3) 

which for genuine, non-radiative DIS events is approximately equal to lE^. Here Ei and j 
are the energy and longitudinal momentum component of a particle i in the hadronic final state. 
This requirement also suppresses events with hard initial state photon radiation. QED Compton 
events are excluded using a topological cut against two back-to-back energy depositions in the 
SpaCal. 

The selection is optimised to obtain large detection efficiency. This required detailed stud- 
ies which were also based on high statistics event samples obtained in the years 2003-2006, 
corresponding to 51 pb~^ of e+p and 45 pb~^ of e~p interactions taken with a dedicated high y 
trigger at 920 GeV proton beam energy. The event selection criteria for the high y region are 
summarised in table 1 . 

The extraction of also requires the measurement of cross sections at lower y. The low y 
region is defined for the 460 and 575 GeV data with y < 0.38 and for the 920 GeV data with 
y < 0.5. The analysis uses a method based on the electron variables for reconstruction and 
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hence is limited to y > 0.1 for all data sets. The data at low y involve large polar angles 6e 

outside the acceptance of the CJC. Therefore in this kinematic region no link to CJC tracks is 
required. At low y the photoproduction background is small and further reduced by a tightened 
cut on Riog < 4 cm. 



2.3 Background Identification and Subtraction 

At low E'^, corresponding to high y, the remaining background contribution after the event 
selection may be of a size comparable to or even exceeding the genuine DIS signal. The method 
of background subtraction relies on the determination of the electric charge of the electron 
candidate from the curvature of the associated track. 

Figure 1 shows the E/p distribution of the scattered electron candidates from e~^p inter- 
actions with the energy E measured in the SpaCal and the momentum p of the linked track 
determined by the CJC. The good momentum resolution leads to a clear distinction between the 
negative and positive charge distributions. The smaller peak corresponds to tracks with negative 
charge and thus represents almost pure background. These tracks are termed wrong sign tracks. 
The higher peak, due to right sign tracks, contains the genuine DIS signal superimposed on the 
remaining positive background. The size of the latter to first approximation equals the wrong 
sign background. The principal method of background subtraction, and thus of measuring the 
DIS cross section up to |/ ~ 0.9, consists of the subtraction of the wrong sign from the right 
sign event distribution in each x, Q"^ interval. 




Figure 1 : Distribution of energy over momentum for tracks linked to clusters in the SpaCal with 
energy from 3.4 to 10 GeV that pass all the cuts listed in table 1. Tracks with a negative charge 
are assigned a negative Ejp. 



The background subtraction based on the charge measurement requires a correction for a 
small but non-negligible charge asymmetry in the negative and positive background samples. 



7 



as has been observed previously by HI [6]. The main cause for this asymmetry lies in the 

enhanced energy deposited by anti-protons compared to protons at low energies. The most pre- 
cise measurement of the background charge asymmetry has been obtained from comparisons of 
samples of negative tracks in e+p scattering with samples of positive tracks in e^p scattering. 
An asymmetry ratio of negative to positive tracks of 1.057 ± 0.006 is measured using the high 
statistics e^p data collected by HI in 2003-2006. This result is verified using photoproduction 
events, with a tagged scattered electron, for which an asymmetry ratio of 1.06 ± 0.01 is mea- 
sured. The difference in the hadronic final state between low and high proton beam energy data 
samples leads to an additional uncertainty of 0.003 on the asymmetry ratio. 

The photoproduction background to the Ep = 920 GeV data, which are analysed at lower y 
than the low Ep data, is subtracted using a PHOJET [15] simulation normalised to the tagged 
photoproduction data. This background estimate agrees well with the corresponding result from 
the wrong sign analysis at high y. 

2.4 Comparison of Data with Simulations 

High statistics Monte Carlo (MC) simulations of DIS events are performed for the three proton 
beam energies using the DJANGO program [16], which includes leading order QED radiative 
corrections. The hadronic final state is simulated using ARIADNE [17], based on the Color 
Dipole Model, with subsequent fragmentation as described in JETSET [18]. The detector re- 
sponse is simulated using a program based on GEANT [19]. The simulated events are subject 
to the same reconstruction and analysis software as the data. The MC simulation uses a QCD 
parameterisation of the structure functions [7] normalised to the measured cross section. 

Figure 2 shows, as an example, comparisons of the 460 GeV high y data with simulated 
distributions, for the energy and the polar angle of the scattered electron prior to and after 
subtraction of the background which is determined using wrong sign data events. The DIS MC 
simulation corresponds to correct sign events with a small contribution from the wrong sign 
events subtracted. The latter are caused by events from lower which can mimic an electron 
cluster at larger and also by charge misidentification for the DIS events at the appropriate 
Q^. The electron energy distribution after background correction is almost uniform. A similarly 
good agreement of the simulation with data has been observed for all other physics and technical 
variable distributions of relevance to this analysis, for all three data sets considered. 

3 Cross Section Measurement 

The scattering cross section is measured in the range 12 < < 90 GeV^ for Bjorken x 
of 0.00024 < X < 0.015. The longitudinal structure function Fl{x,Q'^) is extracted from 
three measurements of cr,. at fixed (x, Q^) but different y = Q'^/sx. The data at lower Ep 
cover the higher y region. In the present analysis the cross section measurement is restricted to 
0.1<y< 0.56 at Ep = 920 GeV and to 0.1 <y < 0.9 at 460 and 575 GeV. 

The measurement of Fl as described below relies on an accurate determination of the vari- 
ation of the cross section for a given x and Q"^ at different beam energies. In order to reduce the 
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Figure 2: Top: comparison of the correct sign data (points) with the sum (open histogram) 
of the DIS MC simulation and background, determined from the wrong sign data (shadowed 
histogram), for the energy E'^ (left) and the polar angle 9e (right) of the scattered electron, for 
the 460 GeV data with E'^ < 10 GeV. Bottom: as top but after background subtraction. 

uncertainty related to the luminosity measurement, which presently is known to 5% for each 
proton beam energy of the 2007 data used here, the three data samples are normalised relatively 
to each other. The renormalisation factors are determined at low y, where the cross section is 
determined by F2{x, Q"^) only, apart from a small correction due to R. Using weighted means 
of cross section ratios, extended over bins at low y, relative normalisation factors are derived 
to be 0.980, 0.995 and 1.010 for the 920, 575 and 460 GeV data, respectively. The relative nor- 
malisation is known to within 1.6%. This uncertainty comprises a systematic error of 1.4%, a 
statistical error of 0.6% and the residual influence of R is estimated to be 0.3%. 

After background subtraction the data are corrected for detector efficiencies and for accep- 
tances using the Monte Carlo simulations. The measured differential cross sections are con- 
sistent with the previous HI measurement [6]. They are shown in figure 3. At large x values 
0",. ~ F2 and the three measurements are in good agreement. The cross sections rise towards 
low X but are observed to flatten and eventually turn over at very low x, corresponding to high 
values of y, where F^ is expected to contribute. This behaviour is consistent with the expecta- 
tion as is illustrated using the cross section as implemented in the Monte Carlo simulation of 
the data. 
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Figure 3: The reduced inclusive DIS cross sections measured at different Q'^ values and shown 
as a function of x for the data taken at the three proton beam energies, 920 GeV (squares), 
575 GeV (stars) and 460 GeV (points). The error bars represent the statistical and systematic er- 
rors added in quadrature. The absolute luminosity uncertainty of the cross section measurement 
is not included in the error bars. Curves for Or as implemented in the Monte Carlo simulation 
of the data are shown as solid (920 GeV), dashed-dotted (575 GeV) and dotted lines (460 GeV) 
while the dashed line represents Fiix^ Q^), which is independent of s. 

The systematic uncertainty on the cross section is derived from various contributions, some 
of which depend on the y region. The uncertainties leading to kinematic correlations are: 

• The uncertainty on the SpaCal electromagnetic energy scale, determined with the double- 
angle method, is 0.4% at large energies degrading to 1% at 3 GeV energy. This is verified 
at the kinematic peak, where E'^ has to be close to E^., and at lower energies with 7r° 77, 
J /^! ^ e+e~ decays and with elastic QED Compton events. 

• The uncertainty on the electron polar angle is 1 mrad, estimated using independent track 
information from the BST, the COZ and the CJC. 
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The hadronic energy scale, calibrated using electron-hadron transverse and longitudinal 
momentum balance, has an uncertainty of 4%. 



• The background charge asymmetry is known to 0.6% based on studies of wrong charge 
data in e^p scattering and tagged photoproduction events. 

• The normalisation of the PHOJET simulation, used for background subtraction in the 
920 GeV data, has a 30% uncertainty. 

• The central track-cluster link efficiency is verified with an independent track reconstruc- 
tion using BST and CJC hit information. The uncertainty of this link efficiency combined 
with the interaction vertex reconstruction efficiency is estimated to be 1.5%. At low y, 
where no track link is required, the remaining uncertainty from the vertex reconstruction 
is 0.5%. 

The uncorrelated systematic uncertainties originate from the Monte Carlo statistical errors and 
from the following sources: 

• The uncertainty on the charge measurement is determined from data to Monte Carlo com- 
parisons at low y and cross checked with radiative events which are background free in the 
low energy region. As the charge misidentification causes signal events to be subtracted 
as background, a 1% uncertainty on Ur is obtained. 

• The radiative corrections are efficiently reduced to below 10 % by the £■ — constraint 
and the topological cut against QED Compton events. A comparison of calculations based 
on the Monte Carlo simulation with the numerical program HECTOR [20] results in an 
uncertainty on cr^ of 1% at high y and 0.5% at low y. 

• The trigger efficiency, determined from independent monitor triggers, is known to within 
1% for the combined CIP-SpaCal trigger and 0.5% for the inclusive SpaCal trigger. 

• Comparisons between different electron identification algorithms and between data and 
simulations yield an estimated uncertainty of 1% (0.5%) on the electron identification at 
high (low) y in the SpaCal calorimeter. 

Further uncertainties, such as the effect of the LAr noise on the cross section, have been in- 
vestigated and are found to be negligible. The subtraction of background using wrong sign 
tracks causes an additional statistical uncertainty which is included in the statistical error. The 
correlated and uncorrelated systematic errors combined with the statistical error lead to an un- 
certainty on the measured cross sections at high y of 3 to 5%, excluding the common luminosity 
error. 

4 Measurement of Fl{x,Q^) 

The longitudinal structure function is extracted from the measurements of the reduced cross 
section as the slope of ar versus y'^/Y+, as can be seen in equation 1. This procedure is illus- 
trated in figure 4. At a given Q'^ value, the lowest x values are generally accessed by combining 
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Figure 4: The reduced inclusive DIS cross section plotted as a function of y'^/Y^ for six values 
of X atQ^ = 25 GeV^, measured for proton beam energies of 920, 575 and 460 GeV. The inner 
error bars denote the statistical error, the full error bars include the systematic errors. The lumi- 
nosity uncertainty is not included in the error bars. For the first three bins in x, corresponding 
to larger y, a straight line fit is shown, the slope of which determines Fl{x, Q"^). 

only the 920 and the 575 GeV data. At larger x, cross section measurements from all three data 
sets are available. These measurements are observed to be consistent with the expected linear 
dependence. 

The central values are determined in straight-line fits to ar{x,Q'^,y) as a function of 
y'^/Y^ using the statistical and uncorrected systematic errors. The systematic errors on Fl take 
the correlations between the measurements into account using an off-set method: all correlated 
error sources, including the uncertainty from the relative normalisation of the cross sections 
which in the extraction of is attributed to the 920 GeV cross sections, are considered sep- 
arately and added in quadrature to obtain the total systematic error due to correlated sources. 
This error is added in quadrature to the statistical and uncorrelated systematic uncertainties to 
obtain the total error on Fl. The measurement is limited to bins where the total error is below 
0.6. 

The measurement of Fl{x, Q'^) is shown in figure 5. The result is consistent with the pre- 
diction obtained with the HI PDF 2000 fit [7], which was performed using only the HI high 
energy cross section data. The measurement is also consistent with previous determinations of 

by HI [6], which used NLO QCD to describe and subtract the F2 term from the measured 
reduced cross section at high y. 

The values on Fl{x, Q'^) resulting from averages over x at fixed Q"^ are presented in figure 6 
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Figure 5: The longitudinal proton structure function Fl{x, Q"^) . The inner error bars denote the 
statistical error, the full error bars include the systematic errors. The luminosity uncertainty is 
not included in the error bars. The curve represents the NLO QCD prediction derived from the 
HI PDF 2000 fit to previous HI data. 



and given in table 2. The average is performed taking the x dependent correlations between the 
systematic errors into account. The measurement of Fl{x, Q'^) is compared with the HI PDF 
2000 fit and with the expectations from global parton distribution fits at higher order perturba- 
tion theory performed by the MSTW [8] and the CTEQ group [9] groups. Within the experi- 
mental uncertainties the data are consistent with these predictions. This consistency underlines 
the applicability of the DGLAP evolution framework of perturbative QCD at low Bjorken x at 
HERA. 



5 Summary 

This letter presents the first measurement of the longitudinal proton structure function in deep 
inelastic scattering at low x. The values are extracted from three sets of cross section mea- 
surements at fixed X and Q^, but different inelasticity y, obtained with three different proton 
beam energies at HERA. The results confirm DGLAP QCD predictions for Q^), deter- 
mined from previous HERA data, which are dominated by a large gluon density at low x. At 
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Figure 6: The longitudinal proton structure function Fl shown as a function of Q'^ at the given 
values of x. The inner error bars denote the statistical error, the full error bars include the 
systematic errors. The luminosity uncertainty is not included in the error bars. The solid curve 
describes the expectation on Fl{x, Q^) from the HI PDF 2000 fit using NLO QCD. The dashed 
(dashed-dotted) curve is the expectation of the MSTW (CTEQ) group using NNLO (NLO) 
QCD. The theory curves connect predictions at the given {x, Q"^) values by linear extrapolation. 



g2/GeV2 


X 


Fl 


Stat. 


uncorr. 


corr. 


total 


12 


0.00028 


0.22 


0.06 


0.05 


0.08 


0.11 


15 


0.00037 


0.08 


0.05 


0.04 


0.09 


0.11 


20 


0.00049 


0.24 


0.04 


0.04 


0.09 


0.10 


25 


0.00062 


0.38 


0.05 


0.05 


0.08 


0.10 


35 


0.00093 


0.24 


0.06 


0.06 


0.09 


0.13 


45 


0.0014 


0.18 


0.08 


0.08 


0.14 


0.18 


60 


0.0022 


0.33 


0.13 


0.13 


0.19 


0.27 


90 


0.0036 


0.48 


0.23 


0.22 


0.22 


0.39 



Table 2: The longitudinal proton structure function Fl{x^ Q"^) measured at the given values of 
and x. The statistical, uncorrected and correlated systematic uncertainties are given as well 
as the total uncertainty. 
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the current level of accuracy, for the covered Q'^ range between 12 and 90 GeV^, the data are 
thus consistent with perturbative QCD. 
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